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We study optical conductivities for high-Tc superconductors under the magnetic field on the basis 
of the microscopic Fermi liquid theory. Current vertex corrections (CVC's) are correctly taken into 
account to satisfy the conservation laws, which has been performed for the first time for optical 
conductivities based on the fluctuation-exchange (FLEX) approximation. We find that the CVC 
emphasizes the oj-dependence of crxy{ui) significantly when the antiferromagnetic (AF) fluctuations 
are strong. By this reason, the relation axyito) ~ {a{Lo)}^, which is satisfied in the extended-Drude 
model given by the relaxation time approximation (RTA), is totally violated for a wide range of 
frequencies. Consequently, the optical Hall coefficient Rii{io) strongly depends on uj below the 
infrared frequencies, which is consistent with experimental observations. We also study the mystery 
about a simple-Drude form of the optical Hall angle Onii^) observed by Drew et al., which is highly 
nontrivial in terms of the RTA since the strong w-dependence of the relaxation time should modify 
the Drude-form. We find that a simple Drude-form of Onii^) is realized because the tj-dependence 
of the CVC almost cancels that of the relaxation time. In conclusion, anomalous optical transport 
phenomena in high-Tc superconductors, which had been frequently assumed as an evidence of the 
breakdown of the Fermi liquid state, are well understood in terms of the nearly AF Fermi liquid 
once the CVC is taken into account. 

PACS numbers: 78.20.Bh, 72.10.-d, 74.72.-h 



I. INTRODUCTION 

In cuprate high-Tc superconductors (HTSC's), vari- 
ous physical quantities in the normal state deviate from 
the conventional Fermi liquid behaviors in usual metals, 
which are called the non-Fermi liquid (NFL) behaviors. 
These NFL behaviors has caused controversial discus- 
sions on its ground state. One of the most predomi- 
nant candidates is the Fermi liquid state with strong an- 
tiferromagnetic (AF) fluctuations. T, ^ 0, 0. In fact, 
spin fluctuation theories like the SCR theory [al and the 
fluctuation-exchange (FLEX) approximation |^ can 
reproduce the Curie- Weiss like behavior of 1/TiT and 
the T-linear resistivity in HTSC's, as well as an appro- 
priate optimum Tc of the order of lOOK with the correct 
symmetry, d^2_y2. 

Especially, anomalous transport phenomena under the 
magnetic field in HTSC's have been long-standing prob- 
lems, as an strong objection against a simple Fermi liquid 
picture. For example, the Hall coefficient Rn is posi- 
tive in hole-doped systems like YBa2Cu307_5 (YBCO) 
and La2-5Sr5Cu04 (LSCO) whereas it is negative in 
Nd2-5Ce5Cu04 (NCCO), although they possess simi- 
lar hole-like Fermi surfaces (FS's) ,7]. In each com- 
pound, i?H oc is observed below Tq ^ TOOK, and 
|-Rh| ^ l/ne [n being the electron filling number) at 
lower temperatures. Moreover, the magnetoresistance 
Po is proportional to T~^ below To 8, 9J. As a re- 
sult, so called the modified Kohler's rule, Ap • po i?^' 
is well satisfied in HTSC's. They cannot be explained 
on the same footing within the relaxation time approx- 
imation (RTA) even if one assume an extreme momen- 
tum and energy dependences of Tk(e): If one assume a 



huge anisotropy of Tk(0) to explain the enhancement of 
i?H experimentally at lower temperatures, then Ap/po 
should increase much faster than experiments (oc T~^) 
because Ap/po is much sensitive to the anisotropy of Tk 
in terms of the RTA. Thus, we cannot explain the modi- 
fied Kohler's rule on the basis of the RTA [Tol[lll|. 

Resent theoretical works have shown that the origin 
of these anomalous DC-transport phenomena in HTSC's 
is the vertex correction for the total current Jk, which is 
known as the back-flow in Landau- Fermi liquid theory ^ 
[T^ . Jk becomes totally different from the quasiparticle 
velocity Vk when strong AF fluctuations exist. Reflecting 
this fact, the DC-conductivities a^i, (p, v = x,y) behaves 
as [12 

a (XT, (T^y oc XQ ■ t'^ , (1) 

where t represents the relaxation time of quasiparticles 
(at the cold spot), and XQ is the staggered susceptibility, 
which follows the Curie- Weiss like behavior in HTSC. As 
a result, i?H oc XQ ^ is concluded. By taking ac- 
count of the back-fiow, we can naturally explain anoma- 
lous behaviors of i?H, the magnetoresistance (Ap/p), the 
thermoelectric power (S) and the Nernst coefficient (i^) 
in a unified way d [II [ll [11 [11 . 

Dynamical transport phenomena in HTSC's are fur- 
thermore mysterious. For example, optical conductivities 
under the magnetic field shows striking deviation from 
the extended-Drude forms in HTSC's. Previous theo- 
retical works, many of them were based on the RTA, 
unable to give comprehensive understanding for them 
,18, 19, 21). In the present paper, we study the 
role of the back-flow in the diagonal optical conductivity 
cr{uj) and the off-diagonal one axyii^) based on the Fermi 
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liquid theory. Here, we develop the method of calculating 
a{uj) and a^yioj) using the FLEX approximation by tak- 
ing the current vertex correction (CVC), which represent 
the back-flow contribution, to satisfy the conservation 
laws. We call it the CVC-FLEX approximation. By this 
approximation, both AC and DC transport phenomena 
in HTSC are explained on the same footing. .221] . 

In the spirit of the RTA, optical conductivities are 
given by the following extended-Drude (ED) forms when 



T ^(w) — iuj ' 



(T-i(tj) - iujy 



(2) 
(3) 



where t{uj) is the relaxation time of quasiparticles. Here, 
w-dependences of and flxy have been dropped for sim- 
plicity. [This simplification will not be allowed in heavy 
fcrmion systems because of the strong w-dependence 
of the renormalization factor.] In usual Fermi liquids, 
T (X {u;^ + (7rT)2)-i. In HTSC, w-dependence of t(w) 
is much stronger; according to a spin fluctuation theory 
r(w) cx {ImTs^^icj -iS)) Fs « (uj + ttT)-^ for a wide 
range of (w, T). Actually, the relaxation time deduced 
from the experimental optical conductivity, which is pro- 
portional to Reo'^^(cj), follows the above relation. 

When the ED-form is satisfied, i?H(w) becomes real 
and w-independent because a^y^{uj) oc {a^^^{uj)}^. 
This relation is approximately satisfied in Cu and Au; for 
w « lOOOcm"^ where w >, t~^{uj) is satisfied, the reduc- 
tion of Rei?H(w) from the DC-value as well as the ratio 
of linRiiii^) to the real part are about 10% for Cu, and 
are about 20% for Au, respectively However, Rni^^) 
shows strong w-dependence in HTSC, which means that 
the ED-forms are totally violated. In fact, we show in 
the present work that axyiuj) strongly deviates from the 
extended-Drude form due to the c<j-dependence of the 
back-flow in the presence of strong AF fluctuations. Ac- 
cording to experiments for the optimally-doped YBCO, 
Imi?H(^) takes the maximum value at wrh ~ 50cm"^, 
and Imi?H(^RH) ~ Rei?H('^RH) ^19.] On the other hand, 
Imcr(tt') takes the maximum value at tUxx ~ 150cm~^. 
This relation wrh ^ i-^xx, which cannot be reproduced 
by the RTA even if one assume strong (k, a;)-dependences 
of Ti^iuj), is well reproduced in the present study. 

We also discuss the optical Hall angle 9b_{uj) = 
axy{oj)/a'{uj), whose ED form is given by 9^"^-^{u}) cx 
{t~^{uj) — iuj)^^. Quite surprisingly, 0-ti{^^) in HTSC's 
follows a simple Drude form even in the infrared (IR) re- 
gion (uj ^ lOOOcm"^) [TslflTl l. For instance, the real part 
of 9^ (lu) in HTSC is almost ^-independent. This exper- 
imental fact cannot be understood in the framework of 
the RTA since the w-dependence of t{uj) is prominent in 
HTSC as mentioned above. Thus, the optical Hall angle 
in HTSC's have put very severe constraints on theories 
in the normal state of HTSC's. In the present paper, we 
show that the simple Drude form of 0h(^) in HTSC is a 



natural consequence of the cancellation between the uj- 
dependence of r and that of the CVC. This fact further 
confirms the importance of the CVC for both DC and 
AC transport phenomena in HTSC's. p^. 

To clarify the reason we derive the general expression 
for the w-linear terms of (t(w) and axy{w) from Kubo for- 
mula based on the microscopic Fermi liquid thoery. They 
are exact up to the most divergent terms with respect to 
T. By analysing the back-flow in the obtained expres- 
sion, we find that the relation lm.axy{uj)/ioj oc ■ r'^ 
holds. Because crxy{0) cx xqt^ cx T~^t'^ according to 
eq.(^, the ED-form in eq.® fails due to the back-flow 
in nearly AF metals. [If we extend DC-axy in eq.Q to 
finite frequencies in the spirit of the RTA, we obtain eq. 
(O with flxy cx xq. However, such an easy extension is 
not true because it gives that l'maxy{uj)/iu> cx XQ'''^-] In 
summary, the enhancement of l'niaxy{i-u)/iLj due to the 
CVC is more prominent than that of axy{0), which leads 
to the violation of the ED-form for f7xy{uj). The present 
study confirms the significant role of the back-flow on the 
optical (as well as DC) conductivities in HTSC's. 

We shortly mention the recent theoretical progress on 
the optical conductivity in Fermi liquids. For exam- 
ple, studies by the dynamical-mean-field-theory (DMFT) 
have been revealed important strong correlation effect on 
the optical conductivity [23. l25| . However, the effect of 
back-flow is totally dropped in DMFT, which is known 
to give the enhancement of i?H in nearly AF metals; see 
eq.(^. We also comment that the effect of back-flow on 
the Drude weight of ^{lj) was studied based on the Fermi 
liquid theory at zero temperature |2^ l^^l . However, the 
overall behavior of 17(0;) at finite temperatures in strongly 
correlated systems is highly unknown. 

The contents of the present paper is the following: In 
§11, we explain how to calculate the self-energy and the 
conductivity by the FLEX approximation In §111, we de- 
rive the exact expression for lim^^—^o 

{Lij)/iuj based on 

the Kubo formula. We discuss the deviation from the 
Fermi liquid like behavior due to the back-flow in the 
presence of the AF fluctuations. In §IV, we address the 
numerical results for <j{u)), axy{oj), Rai^) and 9-a{uj), and 
we compare them with experimental results. We succeed 
in reproducing their characteristic behaviors in a natu- 
ral way at the same time. This is the main part of the 
present study. Summary of the present work is addressed 
in §V. A physical meaning of the back-flow is explained. 



II. NUMERICAL CALCULATION 

A. FLEX approximation for HTSC 

In this subsection, we explain the fluctuation-exchange 
(FLEX) approximation, which is one of a self-consistent 
spin fluctuation theory The FLEX approximation 
is classified as a conserving approximation whose frame- 
work was constructed by Baym and Kadanoff |2^ [2^ . 
In the conserving approximation, correlation functions 
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given by the solution of the Bethe-Salpeter equations 
automatically satisfy the macroscopic conservation laws. 
This is a great advantage of the FLEX approximation in 
studying transport coefficients. In fact, it is well known 
that approximations which violate conservation laws, like 
the relaxation time approximation (RTA), frequently give 
unphysical transport phenomena. 

Origin of anomalous behaviors in the normal state in 
HTSC, which are frequently called the non-Fermi-liquid 
(NFL) like behaviors, have been studied intensively for 
almost 20 years. Recently, many of them are consistently 
explained based on the Fermi liquid picture with strong 
antiferromagnetic (AF) fluctuations, using the FLEX ap- 
proximation, the perturbation theory with respect to [/, 
SCR theory, and so on The range of applicability 

of the FLEX approximation is wide; from the over-doped 
region till the slightly under-doped region {S ^ 10%) 
above the pseudo-gap temperature T* ^ 200K. By tak- 
ing the CVC into account, we can reproduce various 
NFL-like behaviors in transport phenomena above T* 
by the FLEX approximation T^, or even below T* by 
the FLEX-|-T-matrix approximation 14]. As for the or- 
ganic superconductor k-(BEDT-TTF), the d-wave super- 
conductivity 1^ Is^l , as well its Curie- Weiss like be- 
havior of i?H |33 , are also well reproduced by the FLEX 
approximation. 

Here, we study the following square lattice Hubbard 
model: 

H ^k4^Ck<T + »^<T"a : (4) 

kcr i 

where U is the on-site Coulomb interaction, and is the 
dispersion of a free electron. In the tight-binding approx- 
imation, ek = 2t{cos{kx) + cos{ky)) + 4t' cos{kx) cos{ky) + 
2t" {cos{2kx) + cos{2ky)), where t, t' , and t" are the near- 
est, the next nearest, and the third nearest neighbor hop- 
ping integrals, respectively. In the present study, we 
use the following set of parameters: (I) YBCO (hole- 
doped): to = -1, h = 1/6, t2 = -1/5, U = 8. (II) 
NCCO (electron-doped): to = -1, ti = 1/6, = -1/5, 
U = 5.5. (HI) LSCO (hole-doped): to = -1, ti = 1/10, 
t2 = —1/10, U = 5. These parameters are equal to that 
used in ref.|l3|, except that U for LSCO is changed from 
6 to 5. These hopping parameters were determined by 
fitting to the Fermi surface (FS) observed by ARPES or 
obtained by the LDA band calculations. The shape of the 
FS's for YBCO, LSCO, NCCO are shown in ref.'T2]. Be- 
cause to ~ 4000K in real systems, T = 0.01 corresponds 
to - 40K. 

First, we calculate the self-energy numerically using 
the FLEX approximation. The expression for the self- 
energy is given by 

Sk(e„) = Tj2v^'^''''iL0i)Gk-p{en-u;i), (5) 

Lp 

y/LEX(^^) ^ C/'|^X^(c.'0 + ^Xp(^/)-X°(^0|(6) 
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FIG. 1: Schematic k-dependences of Jk and Vk when the 
AF fluctuations with Q ~ (7r,pi) are strong. The cold-spot 
on the Fermi surface for hole-doped (electron-doped) systems 
exist around A (B). 



i-c/x°' i + c/x°' ^ ' 

Xliioi) = -T^Gk(e„)Gk+p(e„+w;), (8) 

where Gk(en) = (*£« + M ^ — Sk(en))^^ the thermal 
Green function, and is the chemical potential. e„ = 
7r(2n -t- l)r, uji = 2'kIT are the Matsubara frequencies 
for fermion and boson, respectively, and x'^ represent 
the spin and charge susceptibilities. We solve the above 
eqs. © and © self-consistently, under the constraint 
of constant electron density n by choosing the chemical 
potential. Hereafter, we study mainly the case for n = 
0.9 for LSCO and YBCO, and for n = 1.15 for NCCO. 
In the present numerical study, 64 x 64 k- meshes and 512 
Matsubara frequencies are used. 

In the present study, the Stoner factor as = C^Xq(O) 
at T = 0.02 exceeds 0.99 both for n = 0.9 (LSCO, 
YBCO) and n = 1.15 (NCCO). (Note that as < 1 at 
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finite temperatures in two dimensional systems because 
Marmin- Wagner theorem is satisfied in the FLEX ap- 
proximation.) In this case, Xq 3> Xq is realized, and 
yFLEX ^ (3C/2/2);^/^ Then, by reflecting the strong 
q-dependence of Xqj 7k = lvaT,i^{—iS) on the FS be- 
comes very anisotropic. (7k represents the quasiparti- 
cle damping rate.) 7k takes a large value around the 
crossing points with the magnetic Brillouin zone (MBZ)- 
boundary, which we call hot spots as often referred to 
in literatures 0, [s^ . On the other hand, 7k takes the 
minimum value at the points where the distance from the 
MBZ-boundary is the largest, which are called cold spots. 
(see Fig. ^) Transport phenomena for lower frequencies 
are governed by the electronic properties around the cold 
spot. 

As shown in figQ] (a), the location of the cold spot for 
hole-doped systems is around A. Whereas the cold spot 
for electron-doped systems locates around B: This fact 
was first predicted by ref. ^3 theoretically based on the 
FLEX approximation, and it is verified by the ARPES 
measurement later jSy, [s^l • The difference of the posi- 
tion of the cold spot explains the remarkable difference 
in transport phenomena between hole-doped systems and 
electron-doped ones T^]. For example, both i?H and S 
are positive in YBCO and LSCO whereas they are nega- 
tive in NCCO. 

In the present FLEX calculation for LSCO {n = 0.90), 
the maximum (minimum) value of 7k on the FS is 0.38 
(0.12) at T = 0.02. The ratio of anisotropy is 3.2, which is 
consistent with ARPES measurements in slightly under- 
doped compounds. Such a small anisotropy cannot ac- 
count for the enhancement of i?H in under-doped systems 
within the RTA ■ In the presence of strong AF fluc- 
tuations, however, the total current Jk becomes quite 
different from the quasiparticle velocity Vk due to the 



back- flow, which is totally dropped in the RTA l^. In 
fact, various anomalous transport phenomena in HTSC's 
are brought by the nontrivial momentum-dependence of 
Jk around the cold spot. 

Finally, we discuss the dynamical spin susceptibility 
given by the FLEX approximation. Its phenomenological 
form is expressed as 



XQ 



i + e'(q-Q)2-«^M: 



(9) 



where ^ is the AF correlation length, and Q is the nesting 
vector. Q = (tt, tt) in both YBCO and NCCO JJj. XQ « 
l/cjgf oc 1^^, and oc in the FLEX approximation, 
which is equivalent to the SCR theory 0. Here we call 
eq.© the AF-fiuctuation model. 



B. Conserving approximation for (T^^(a;) 

Here, we explain how to calculate the optical conduc- 
tivity based on the Kubo formula. To satisfy the conser- 
vation laws, we have to take the current vertex correction 
(CVC) into account in accordance with the Ward iden- 
tity. 

According to the Kubo formula, the optical conduc- 
tivity axxi^) = f (t^) and the optical Hall conductivity 
axyiijj) are given by 



1 



(10) 



where K^^{uj) is the retarded current-current correlation 
function, which are given by the analytic continuations 
of the following thermal Green functions Issl l3fll | : 



Kxx{ii^i) = -2e^r^w£^gk(e«;'^i)Akx(en 

KxyiiUJl) 



(11) 



n,k 



n,k 

If d 

X- < [Gk{en +i^i)vkx{en +1^1) - Gk{en)vkxien)]-Q^ (x ^ y) ^ Aka(e„;w;), 



(12) 



where — e(e > 0) represents the charge of an electron. 
Here we put c^l. gk{en;uji) = G'k(en -t- w/)Gk(e„), 
v^^j = de^^/dk^ is the velocity of the non- interacting elec- 
tron with k, and Ak;^(eri; ^i) is the dressed current which 
is given by 



Ak(e„;w;) = + r^rk,p(e„,e„i;a;;)gp(e„i;a;;)v° 

m,p 



m,p 



xAp(e„; Wi), 



(13) 



where rk.p(e„, e,„; wj) is the full four-point vertex func- 
tion, and rjj. p is the irreducible four-point vertex with 
respect to the particle-hole channel. 

Equation l(TT|l was derived by Ehashbcrg in 1962 [38l |. 
Equation 112|) was first derived within the Born approx- 
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imation |40j. later it was proved to be exact for any 
Fermi liquid system up to the order of 0(7"^), which is 
the most divergent term with respect to 7^^ 39]. In the 
same way, the general formulae for the magnetoresistance 
the thermoelectric power and the Nernst coefficient 
j42| have been derived from the linear-response theory 
(Kubo formula) on the basis of the microscopic Fermi 



liquid theory. 

In the numerical study of transport coefficients, one 
have to take account of the irreducible four point ver- 
tex to satisfy the conservation laws, which is given 
by 6J:^^^^/6G 28]. In the FLEX approximation, F^ is 
expressed as 



-T^ (G(fc + k'+p) + G{k - p)) Gi-p + k')Wip;p + coi), (14) 

pm 

W{p;p + oji) = ^U^{l + Ux%p))il + UxHp + iOi)) + ^U^{l-Ux'{p))il + Ux%P + ^i))-U^ (15) 



where k,k' ,p = (k, e„), (k', e„'), (p, a;;). In the light hand 
side of ea. (|14|l . the first and the second terms are called 
the Maki-Thompson (MT) and the Aslamazov-Larkin 
(AL) terms, respectively, in literature. 

Using the Green function given by the FLEX approx- 
imation, we numerically solve the Bethe-Salpeter equa- 
tion for A, ca. H13l) . by iteration. The kernel function F^ 
is given in ea. (|14|l . Then, we obtain K^^{iuji) by insert- 
ing the obtained A into eas. Hll|) and H12|) . The retarded 
function K^^{u}) is derived from the analytic continua- 
tion of Kf^^{iu!i) with LUi > 0, using the numerical Fade 
approximation. 

Fade approximation is less reliable when the function 
under consideration is strongly w-dependent, and when 
the temperature is high because the Matsubara frequency 
is sparse. To increase the accuracy of the Fade ap- 
proximation, both S(e„) and A(e„;e/) have to be ob- 
tained with high accuracy; their relative errors should be 
10~^ ~ 10~^°. In performing the Fade approximation, 
we utilize the fact that the ia;-lincar term of Kf^^^uj) is 
equal to the DC value of a^,^, which can be obtained 
within the FLEX approximation with high accuracy, as 
performed in refs. Ufll, 12, 13, 14]. By imposing this 
constraint on the Fade approximation, we succeed in de- 
riving the a^^iuj) with enough accuracy in the present 
study. 

In the present numerical study, we take the infinite se- 
ries of the MT-terms in A, whereas we drop all the AL- 
terms in ea. (|14|) . This simplification is justified for DC- 
conductivities when the AF fiuctuations for Q « (tt, tt) 
are dominant, as proved in ref. 12]. In the same way, 
AL-terms would be also negligible for IR optical conduc- 
tivities. In fact, we have checked that the /-sum rule 
both for a{uj) and axy{Lo) are well satisfied even if all 
the AL-terms are dropped, as will be shown in FigllOl 
This results ensure the reliability of the present numer- 
ical study. (Although we have also tried to include the 
AL-terms, then the accuracy of the numerical Fade ap- 
proximation became worse, unfortunately.) 



In a Fermi liquid, the relaxation time Tk in the 
relaxation time approximation (RTA) is l/27k = 
l/2ImEk(— 1(5). Transport coefficients can be expanded 
in terms of 7~^, which diverges as the temperature ap- 
proaches zero. Equation Hill) is formally an exact expres- 
sion. On the other hand, ea. (|12|l gives an exact expres- 
sion for (Txy{uj) up to the order of 0(7^^), which is the 
most divergent term with respect to 7^^. Less singular 
terms, which are given as C and D in ref. ,32] (p. 636) or 
the last two terms in Fig. 3 of ref.^^li ^i'^ dropped in 
Kxy given by eq. H12|) . In terms of the FLEX approxima- 
tion, ea. ((T^ is exact beyond 0(7"^) because both C and 
D vanish within this approximation js^l- 



III. c^-LINEAR TERM OF (7^^(cj): THE ROLE OF 
THE CVC 

Based on the Kubo formula, we derive the general 
formula for the w-linear term of a^i,{u)), which we de- 
note Acr^^ hereafter. The relation 17^^(0;) = {a^^^—Lo)}* 
tells that Acr^jy is pure imaginary. The most divergent 
terms of Act = t^a^x and Aa^^y are order of 0(7"^) and 
0(7^^), respectively. The derived expressions in this sec- 
tion are exact up to the most divergent terms. Based on 
the derived expression, we discuss the role of the current 
vertex correction (CVC) for Act^j^, and find that it is 
strongly enhanced when the AF fluctuations ate strong. 
Readers who are not interested in the microscopic deriva- 
tion can skip to the next section, where we will show the 
numerical results for af^^{uj) given by the CVC-FLEX 
approximation. 



A. Exact expression for the oj-linear term of an^{(jj) 

Here, we perform the analytic continuations of eq. Hll|l 
and (|12|l according to refs.js^ and [32|], and derive the 
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general formula for AtT^^(w). In preparation for analyz- 
ing the CVC in later sections, we derive the expressions 
without CVC, which corresponds to the relaxation time 
approximation (RTA) with k- and e-dependent relaxation 
time Tk(e). After the analytic continuation of ea. (|ll|l . 



de 1 

TT 2UJ 



th 



"2T 



th- 



2T 



xG^(e+)G^(6-){<}^ 



where e± = e±a;/2, and we have dropped the terms with 
G"^(e+)G^(e_) and G^{e+)G^{e-) because they are less 
divergent with respect to 7~^, and they vanish at T = 
[3^. Here we expand ea. (fTC|l with respect to uj: 



L--th- = 2c. 

R/, ^/^A/ \ _ 1^ / \|2 



(17) 
(18) 



x|Gk(e)|2|ImGk(e)|+0(^' 



where z\f^[e) = [1 — ^ReSk(e))^^ is the renormalization 
factor. At sufficiently lower temperature, the Green func- 
tion for Lj ^ Q and et ^ is well approximated as 



1% 



(19) 



where ej^ — 2:k(ek + R.el]k(0) — /i) and — Zk7k- Equa- 
tion is called the quasiparticle representation of Green 
function, whose validity is assured by the microscopic 
Fermi liquid theory. When ea. (|19|l is valid, 

|ImGk(a;)| = 7rzk(5(tJ - e*^), 



|Gk(^)|' 



7k 



Gk(c^)P|ImGk(6)| = -^z^5{uj - ej;), (20) 



27, 



|Gk(^)|^ = — Zk<5(c.-4), 



TT 



|Gk(^)|^ReG^(e) 



'47^ 



Zk<5(w - 



for (jJ^e'^ ~ 0. As a result, the expression for o""(a;) ex- 
panded with respect to w is 



de 



7k 



(21) 



dej 27^ 

In a free-dispersion model Ck ~ k"^ /2ra, ea. (|21|l becomes 
a«(c.) = 



m • 27 



+ iz ^LO 



TO • (27)2 



2 

en 



' to(27-zz-M' ^^^^ 
which is equal to the Drude form for cr(w) given by the 
RTA if we replace 27 with r~^. 

In the same way, <Jxy{'jj) without CVC is given by the 
analytic continuation of ea. lfl^ : 



(16) ^ly{^) 



TT 2uj 



'2T' 



x[G^(e+)-G^(6_)] -Ak, 







(23) 
(24) 



where we have dropped CVC for the quasiparticle veloc- 
ity given by the momentum derivative of the self-energy. 
Up to the order of 0{uj), we see that 



GlGt [G« -G^]^ = |Gk(e)nimGk(6) 



z^'u:\-\G^{er^\G^{erReGi{e] 



— 2 ^5{e ~ el) -f iz^ ^^Trr^Sie ~ e*^). 
^7k ^7k 



(25) 



As a result, the expression for (j'^y{i-o) within the order of 
0{uo) is 



3^ / df\ Ak 



de)2^l 



de 



27^ 



In a free-dispersion model, ea. H26l) becomes 



+ 2iz~^uj- 



m? ■ (27) 



TO^ • (27)'^ 



m?{2^ — iz ^LoY ' 



(27) 



which is also equal to the Drude formula given by the 
RTA. 

In the next stage, we derive the general expression for 
^(^tiu by taking all the CVC's into account. After the 
analytic continuations of eas. (|ll|l and (|12|l [sM l39|. 
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cr(w) 



(28) 



r 



{x y) 



>/ky(e;tj), 



(29) 



where g^^\e;u;) = G^(e+)G« (e_), g['^(e;^) = 
G«(e+)G^(e-) and ^(^^(e;^;) = G^(e+)G^(e_). 
y^^(e) = <^ + ^^^{e). and G'^ are the retarded 



and advanced Green functions, respectively. J, v and v 
are given by 



Jk(e;w) = Vk(e;t^) + ^ J —T22{-ke,(ie';Lu)g^^\e';uj)vi,{e';Lu) 

= Vk(e;^) + X^ J ^4°)(ke,qe';^)5(2)(e';a;)Jk(e';u;), (30) 
Vk(6;^) = vO+ 5] /^r^f (ke,qe';c.)5«(e';c^)v0, (31) 

Vk(6;.) = E -2)/ ■ th + ^''^y^^"^^ ) 4^)(e-;..)r,-.(q6Mc.;^), (32) 

I 



where the definition of 7^2°'' (^e, qe'; ^) is given in ref . [38l| . 
7^^2 is a subgroup of which is irreducible with respect 

(2) (1 3) 

to (7q , whereas it is reducible with respect to gq' ■ The 
following Bethe-Salpeter equation holds; 722 = T^^i + 



E/^2?^5(^^r22. 

As explained in ref.js^, Wk2;(e;w = 0) = v\^x{^\^ — 
0) = Reiiy (£)_ = v \ {^) is well satisfied in a Fermi liquid. 
As a result 0,^, 



a(0) 

CT:ry(0) 



^'E/vf-f ) |Gk(e)|V.(6)Jk.(e), 



TT 



(33) 



-^'E / V (-^) |Gk(e)nimGk(e)||vk(6)|Jk.(e)^Jk,(6) 

k ^ II 

9Jk(e) 



-^'E / ^ (-^ ) |Gk(e)P|IniGk(e)||vk(e)| ( Jk(e) x 



(34) 



where fcy = (e^ x v)/|v|, which is parallel to the Fermi 
surface. In deriving the first line in ea. H34() . we have used 
the relation gf- - Wk^ gf^ = (e. x v)V = |v|gfj. The 
Onsager's relation o^y ~ —Oyx is used in deriving the 
second line in eg . (|34|l . 



Here we expand Jk(e;w) with respect to uj as Jk(e = 
Q;uj) = Jk(0) + icjJ^^^(O). From eq.^, one can check 
that the most divergent term of i''^^ is proportional 
to 7~^, which comes from the w-derivative of g'^^-' or 
that of the thermal factor in 7^^°''; ^^2°^ (^i ^'s ~ 



8 



2(-a//9e')Rer(0,0) + 0{lu^). For simplicity, we de- 
note hereafter J{e;uj = 0) = J(e), r(ke, pe'; cj = 0) = 
r(ke,pe'), and so on. On the other hand, the oj-hnear 
term of ea. H^ll) or (|32|l is not singular with respect to 
7"-^, so we put = in eas. (|^ and hereafter. 
Using the relations in ea. (|20|l . jj^^' oc 7^^ is given by 



r(i) 



y- 

q 



-r,,(ko,qe;)^.^ 



1^ 

7q 



Jc 

27^ 



x5:Rer(q,q') i-f^ 

q' 

Lk + Mk, 



27q' 



(35) 



where the first and the second terms come from the lo- 
derivatives of g^^^ and 'T!^\ respectively. Then, the w- 
linear term of the conductivity, which we denote as Acr^j/ , 
is given by 



Bk = 



27^ 



7k 





■^k^^k Sk" 




[ H 27^J 



JkxIVkl 



dk 



-Jky 



d 

4LVki^Jk, 



(1) 



Jk.|vk|^Jk, 



(36) 
(37) 

(38) 
(39) 



where |vk|gfj = Wkxgf^ - "kygf;^ as explained before. 
They are exact with respect to the most divergent terms 
with respect to 7~^. Note that Acr^,y is pure imaginary. 

Here, we further analyze 3^^^ given in eq. (|35|l . First, 
Lk is given by the following Bethe-Salpeter equation: 



27i: 
Jk 

27k 



E^^22(k0,qe* 



q 

:^l+Ei^22(kO,qe;)^.P„ 



7q 27* 



q 



27k 



7q 



(40) 
(41) 



where TiMqe*^) = 2iImV«„q(e* )(cth(e* /2T) - 
th(e*/2T)) in the FLEX approximation. In a similar 
way, Mk is rewritten as 



1 



^ 4i 
q 



7q 



Nk^-El^er(k,q)(-|)^^.,A. 



(42) 
(43) 



At sufficiently lower temperatures, the expression for 
A(T^^(w) in ea. H37|l is rewritten by using J^^-* (instead of 
jW) as. 



Act = 



7k 



k 



dl 
de 



27^ - 



"^k 



4LVki 



dk 



9 



(1) 



Mk = Jk 



27k' 



(44) 

(45) 
(46) 
(47) 



In the next subsection, we will discuss the temperature 
dependence of /S.axy when the AF fluctuations are strong, 

by analyzing the k-dependence of Jk^^ We will approxi- 
mately solve the Bethe-Salpeter equations (|40|l and H42(l 
based on the AF-fluctuation model. 

In a Fermi liquid, 7k and Zk are expressed as 

^ ^ ^der2i2(kO,qe)pq(e) 



7k 



q 



lE^22(kO,q<). 



(48) 



d 



ReSk(e) 



de 

ERer(k,q) 
q 

-AReEk(O). 



e=0 



5/ 
de 



(49) 



Note that the uniform charge susceptibility in a Fermi 
liquid is given by Xc = [1 — c?ReSk(0)/9^]x°, where Xc 
is the susceptibility for U — Q. Because Xc ^ Xc is 
expected in strongly correlated systems (like in heavy 
Fermion systems), the relation 9ReI]k(0)/9/i <, 1 should 
be satisfied. This relation is also expected to be realized 
in HTSC according to the FLEX approximation. 



B. Role of the CVC in the presence of AF 
fluctuations 



Using the general expression for Acr^,y derived in the 
previous subsection, we discuss its temperature depen- 
dence when the AF fluctuations are strong. For that 
purpose, we approximately analyze the CVC included in 
the expression for cr^;^(a;) based on the AF fluctuation 
model given in eq.©. 

First, we explain the total current for the DC- 
conductivity. The Bethe-Salpeter equation ea. l|^ is 
rewritten at lower temperatures as |l2l| 



Jk = Vk + E^22(kO,pe;)zp- — -Jp (50) 

7p 



p 
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for e = Lu = 0. In the FLEX approximation, 
T22(k0,pe) « 3t/2.iInix^k_p(e+i(5)[cth(e/2T)-th(e/2T)]. 
Due to the thermal factor, 722(kO, pe) takes large value 
only when |e| <, T. If we apply the AF-fluctuation model, 
eq.©, the main contributions of the p-summation in 
eq.(jnni) come from the region |p — k'| <, 1/^, where 
k' is the momentum on the FS defined as {k'^,k'y) — 
—sgn{kxky) ■ {ky, kx), as shown in Fig. ^(b)- We see the 
relation k — k' « Q is satisfied on the FS. Here we as- 
sume IQ — (k — k')| ;< even at the cold spot, which is 
in fact satisfied in the present FLEX approximation for 
hole-doped systems In this case, 7k^ cx T is satisfied. 

Taking account of the expression for 7k given in 
ea. (l48ll . we obtain a simplified Bethe-Salpeter equation 

m, 



In fact, in the present FLEX calculation for LSCO (n = 
0.90), the maximum (minimum) value of ImSk(— «(5) on 
the FS is 0.38 (0.12) at T = 0.02; the ratio of anisotropy 
is 3.2, reflecting the sharp q-depcndencc of lmXq{uj) / uj . 
On the other hand, the maximum (minimum) value of 
z^^ — 1 on the FS is 5.0 (3.8); the ratio of anisotropy is 
only 1.4. 

Considering ea. (|49|l . we rewrite ea. (|55|l as 



Jk' _ / -1 



Xc/x°) ^ -5k;-^. (56) 

^7k 



<5k ;S ttk ;S ttk should be satisfied because Xc/Xc P'^S" 
itive. We note again that Xc/Xc ^ 1 in strongly corre- 
lated systems. Then, the approximate solution of eg . H42|l 



Jk = Vk + ak • J 



k'- 



(51) 



where ak = \Y.p'^22i'^^^P'^p)zp^p/lp\/\^i<i\- Accord- 
ing to the AF-fluctuation model, ak ~ ( cos(0j(q) — 
^j(k')) )|q-k'|<i/c « (1 - c/e) < 1 where c ~ 0(1) 
is a constant, ak takes the maximum value around hot 
spots. The solution of eq. (|51|l is jl^ 



1 



1 



(vk + ak ■ Vk') , 



(52) 



whose schematic behavior is shown in fig. ^ (b). 
Note that {vk'x,Vk'y) = -sgii{kxky) ■ {vky,Vkx) and 

(Jk'x, Jk'y) = -Sgn{kxky) ■ (Jky, Jka;)- 

We stress that the same vertex functions 722(kO, pe) 
appear in eas. H30|l and (|48|) . which is the consequence 
of the Ward identity and is satisfied in any conserving 
approximation. This fact assures that the relation 1— a cx 
f ^ holds even beyond the FLEX approximation. 



In the same way, we study J 



Pk + Mk. The 



Bethe-Salpeter equation for Pk, eq. (|40|l . is simplified as 



Pk-^ + akPk' 
27k 



(53) 



where ak is the same as that in ea. (l51|l . The solution is 
given by 



1 Jk + QkJk' 

27^ l-al 



(54) 



Next, we analyze Nk in ea. (|43|l . Considering the rela- 
tion ReFk.q ^ z-^ReV^}:^^{0), we obtain 

Nk « -^•ak^Rer(k,q) (^-^^ Zq, (55) 



where ak ^ ak is expected in general because the mo- 
mentum dependence of ImV^j'"p''^(e)/e is much promi- 
nent than that of ReVifi^^^(O), which are included in 
eas. (|^ and (j^SJ, respectively. Actually, ImXq(i^)/a; cx 
{Xq(0)}2 cx ^4(1 -t- ^2(Q _ q)2)-2 according to eq. ®. 



Mk 



iSk QkJk + Jk' 

"27i: 



(57) 



In conclusion, an approximate expression for axy{i^) 
up to the order of 0(a;) is given by 



Ca;j/(w) 



X Jk(0;w) X 



dkn 



(58) 



Jk(0;cj) = 3k{0: iv) + luj-^ =J^ + iLU j'^K (59) 

27k 

where j'^^ is given by eqs. l|T7|) . ifS^ and (|S7|l . After a 
simple but lengthy calculation, J x i9J/i9fc|| in ea. (|^ is 
rewritten as 

^ /n ^ 9Jk(0;tj) aJk 
Jk(0;w) X 5^ = Jk X 



dki 



iu! 1 — akO^k 



Jk X 



9/fc|l 

5Jk 

dk\i 



7k l-«k 

9/lak-ak\, 

+^^^ TT^-, — Jk X Jk' 

5^11 V27k l-"k / 
+0{u;'), (60) 

where the second and the third terms contribute to Aa^y 
Note that 



Jk X 



dkn 



,2 1 del 



9 9 

(Jk X Jk')^ = sgn(fc^fcy)-Y^— (62) 

where 9'^ — tan^^{Jy/ Jx)- At the cold spot in hole- 
doped systems, ea. l|^ is proportional to cx be- 
cause |Jk| <, |vk| and [dO'' /dk\\) cx • {d9^/dk\\) around 
the cold-spot in hole-doped systems (point A in Fig^l, 
which we denote as kc hereafter. Note that [dO"" /dk\\) 
represents the curvature of the FS p^ . 



10 



Let us consider the hole-doped system, where the cold 
spot locates on the XY line in Fign(a). Because Jk x 
Jk' = at the cold spot, the second term of the right- 
hand-side of ea. (|60|l gives the main contribution to Acr^^j,. 
This fact immediately tells that 



(63) 



whereas a = = 1 if all the CVC's are dropped, i.e., in 
the RTA. kc represents the cold spot. Considering ea. l|Bn|) 
and the relation 1 — a oc we obtain that m « 4 when 
ttk ^ ttk, and m « 2 when fik ~ ctk- As discussed above, 
fSk ^ ctk is expected by the present analysis for the CVC. 
In the next section, we will show that the relation m w 3 
holds for hole-doped systems in the numerical study. 

In a similar way, we also discuss the role of the CVC 
in a{Lu) for hole-doped systems. Because Jk — — Jk' at 
the cold spot kc, we obtain 



1 + 5k, Jke 

1 + ak, 27*^ ' 



(64) 



which is close to Jko/27k, because ak, , Ok, £ 1. We note 
that |Jk, I < |vk„| because Jk„ « Vk,/(1 + a^J. This 
result suggests that the CVC changes the values of a 
and Act only slightly. As a result. 



cr(aj) cx Wko-^kc 



1 



27k. 



iz ^oj 



1 



(27k J 



(65) 



In conclusion, (t(w) is insensitive against the CVC, as 
the case of DC-conductivity within the FLEX approx- 
imation We will show in the next section that 
cr(w) ~ a^'^^{uj) holds for hole-doped systems in the nu- 
merical study. 

According to eqs. H63() and (|65|l . the Hall coefhcient 
and the Hall angle are given by 



-Rh(^) oc a + iujz ^ {b 

Pjj \UJ) cx ILOZ 



1 



a 



7k, 
26 1 

/^2 r. 



(66) 
(67) 



where a = 5 = 1 in the absence of the CVC, that is, in 
the RTA. 

As a result, we can conclude that the origin of the 
anomalous behaviors of i?H(^) and dYiii^)^ that is, promi- 
nent deviations form the extended-Drude (ED) formula, 
is the strong temperature dependences of a and h which 
originate from ^. We will discuss this mechanism in more 
detail in later sections. 



IV. NUMERICAL RESULTS 

In this section, we show the optical conductivities ob- 
tained by the FLEX approximation with full MT-type 



2 



LSCO (n=0.9) 



full CVC's 
CVC by dE,/dk, 
without any CVC's 



0.05 



0.1 

T 



0.15 



0.2 



LSCO (n=0.9) 

full CVC's 
CVC by clE,/dk„ 
without any CVC's 




0.05 



0.1 

T 



0.15 



0.2 



FIG. 2: Resistivity and the Hall coefRcient for ctJ = ob- 
tained by the FLEX approximation: (i) without any CVC's, 
(ii) with CVC given by dEk(0)/dfc^, (iii) with all the CVC's 



(CVC-FLEX). At lower temperatures, i?H 
and i?H oc T^^'"^ for (i), respectively. 



for (iii) 



CVC's. This kind of calculation has been performed for 
the first time. Hereafter, we call this scheme "the CVC- 
FLEX approximation". We will see that Oxyiy^) shows 
striking deviation from the ED-form, which is highly con- 
sistent with experimental results. Here, the unit of en- 
ergy is the nearest- neighbor hopping integral ^Oi which 
corresponds to ^ 4000K according to LDA band calcula- 



tion. Thus, T ^ 0.01 
the present study. 



40K and w = 0.1 - 300cm- 



A. DC transport coefficients 

Before discussing the optical conductivities, we shortly 
explain the DC transport phenomena given by the CVC- 
FLEX approximation . Obtained p and i?H are shown 
in Fig|21 Results by the CVC-FLEX approximation, 
which are calculated using eas. (|33|l and (|34|) . are denoted 
as "full CVC's" in figures. Curie- Weiss like behavior of 
i?H (more precisely i?H oc T~"-®^) is reproduced due to 
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the CVC. In the Fermi liquid theory, the CVC is divided 
into (i) the back-flow which is expressed by 7^2 j and (ii) 
the renormahzation of Vk given by dTj\^/dkx- To clarify 
the effect of the back-flow, we calculate the conductivi- 
ties by replacing all the Jkp(e)'s with Uk;i(e)'s in eas. l|5^ 
and The obtained results are denoted as "CVC by 

dEk/dfcj;" in fig 13 They correspond to the "without VC" 
in ref. J^]. We see that the resistivity increases to some 
extent due to the back-flow (722)- 

Furthermore, we calculate the conductivities by replac- 
ing all the Jk^(e)'s and WkAt(e)'s with in eas. (|S^ and 
(|34|l . The results are shown as "without any CVC's" in 
fig[21 Then, the resistivity takes the smallest value be- 
cause self-energy correction for the velocity enhances the 
conductivity; |vk(0;0)| > |vj^| at the cold spot. Here- 
after, "RTA" in figures represents the results by "with- 
out any CVC's". As will be shown in figs. 2|and|31 the 
DC-conductivity by RTA is smaller than that by the con- 
serving approximation, because the effect of the velocity 
correction dominates the back-flow effect. 



B. and axy{Lo) 

Here, we perform the numerical calculation for the 
complex optical conductivities, (t{uj) and Oxyioj), using 
the FLEX approximation. The CVC is taken into ac- 
count in the conserving way. We calculate (Jfiu{ijo) by 
ea. (|10|l . where K^^{lo) is derived from eqs. and H12|) 
using the Fade approximation. As explained above, we 
utilize the values oi ij{uj = 0) and Uxyi'^ = 0), which are 
derived from eqs. (|33|l and H34(l as shown in ref.^3|j in 
the course of the Fade approximation. This procedure 
is highly demanded to achieve enough accuracy. 64 x 64 
k- meshes and 512 Matsubara frequencies are used in the 
present FLEX approximation. 

Here we derive an extended-Drude (ED) forms for 
CT°j,(w) from the Kubo formula within the RTA, where 
the suffix means the result by the RTA hereafter. At 
zero temperature, <y'^{uj) for smaller uj (i^-' <, 7) is given 
by, 



e 

UJ 



E 



Z o 



iV(0)- 



i7k(e) 



rfe- 



2< 



7kc (e) + 7kc (t^ + e) - i^k, ^ 



(68) 



where k^. represents the cold spot, ijj-dependence of z{lo) 
has been neglected. In deriving ea. H68|l . we take only the 
contribution comes from the cold spot into account. In 
the same way, 



-iV(0)- 



real 



LSCO (n=0.9) 




0.06 

0.1 

RTA(T=0.02) 

3 O Sy^pCo)) 




LSCO (n=0.9) 



real 



— u 
O 





imaginary 



T=0.02 
0.04 

0.06 

0.1 

RTA (T=0.02) 
6Yed(™) 



0.2 



0.4 



CO 



FIG. 3: Obtained w-dependence of (a) {(t(lj)} ^ and (b) 
{a^y(uj)Y'^'^ by CYC-FLEX approximation. Suffix '0' repre- 
sents the results by RTA. Real parts of {cr(aj)}"\ 
and {a\y{uS)\~'^'^ are approximately proportional to 7ed(i^) 
for Lo < 0.2, as expected by the extended-Drude expression. 
However, tj-dependence of Ke{axy{Lo)}~'^'^ is much weaker 
than that of 7ed(<j-') due to the CVC. 



de- 



(7k, (e) + 7k, (w + e) - i^k,^w) 

(69) 

where Ak is given in ea. (|24|l . In a crude expectation, 
7kc(e) + 7ko(w -I- e) in eas. (|^ and would take the 
minimum value around e ~ —uj/2. As a result, we ob- 
tained the following ED expressions for smaller uj: 



ED/ \ 

a [uj) 



n 



27ed(w) -I- ^uj' 



n 



" (27ED(t^) + «z 



+ IZ ^UJ) 

where 7ED(t^)(> 0) is approximately given by 
1 



7ed(w) 



[7k,(u;/2)-H7k,(-c^/2)]. 



(70) 
(71) 

(72) 



Below, we will show that the above ED formulae for <j{lli) 
still holds even if CVC is taken into account, whereas it 
completely fails for (Jxy{i^) owing to the CVC, which is 
the origin of anomalous behaviors of R}i(lo) and Oh^lo). 

Figure |31(a) shows the obtained cr-^u) for LSCO. We 
see that Re{a~^{uj)} possesses strong w-dependence so 
the simple Drude formula is violated. It is approximately 
proportional to 7ed(^) for lo <, 0.2, which suggests that 
the ED-form in ea. (|70|l is well satisfied, even if CVC is 
taken into account. We also see that lm{a{uj)}~^ shows 
moderate uj- and temperature-dependences, which is pro- 
portional to u!z~^{uj) according to the ED-model. As 
shown in fig. O its gradient decreases as w and/or T 
increases, which is naturally explained as the lo- and T- 
dependences of z~^{uj) (sa m* /ra). 

Figure 121(b) shows {a^yiuj)}-^"^ for LSCO. We rec- 
ognize that Re{(T°j^(a;)}~°-^ cx 7ed(w) within the RTA, 
whereas Ke{axy{oj)}~'^'^ given by the CVC-FLEX ap- 
proximation possesses much moderate tj-dependence. 
These results means that the ED-form in eq.ljjlj is well 
satisfied for a^y{uj), while it is violated for <Txy{u!) due to 
the CVC. We will study the role of the CVC in (Txy{uj) in 
more detail hereafter. We also see that Imjcr^y (tj)}^"-^ 
is almost unchanged against the temperature, whereas 
its gradient slightly decreases as uj increases. 

Figures 01 and [S] show the cj-dependence of <t{uj) for 
LSCO {n = 0.90) and YBCO {n = 0.90), respectively. 
Both of them are qualitatively similar. Parameters for 
each compound are explained below eq.Q). In both 
cases, a{uj) given by the CVC-FLEX approximation is 
slightly larger than that by RTA. In more detail, 
decreases due to the back-flow, whereas it increases due 
to 9Ek/9fc^ in Vk(e); the latter shghtly dominates in 
the present model parameters. Rea{u}) apparently de- 
creases much slower than Lorentzian for larger u, be- 
cause Jed{'-^) increases with uj. The Drude weight is a 
little sharper in LSCO because a smaller value of U is 
used. We note that the Drude weight increases as the 
system moves away from the half- filling (n = 1). For 
LSCO at T = 0.02, Ima(a;) takes the maximum value 
at iOxx ~ 0.06, which is about three times larger than 
27kc(0) because of the w-dependence of 7k (w). 

In contrast to <j{llj), axy{uj) with full CVC is quite 
different from cr°j^((jj) given by the RTA: Figures El 
and 13 show the w-dependence of axy{u>) for LSCO and 
YBCO, respectively. The w-dependence of Recr:Ey(^) be- 
comes prominent due to the CVC. For LSCO (YBCO), 
Keaxyit-o) at T = 0.02 takes a large negative value for w > 
0.03 - 120cm-i {oj > 0.045 ~ ISOcm"^ which is consis- 
tent with experimental observations [isl IT^ UtI IT^ . 
Although Re(T°y{uj) also changes its sigh for ui > 0.1, 
its absolute value is very small. It is naturally under- 
stood from the ED-form because 7(w) increases with uj. 
This large dip in Iieaxy{uj) is naturally understood in 
terms of the /-sum rule, ea. lj55|) . because Reaxy{0) > 
takes an enhanced value due to the CVC. We also stress 
that l'niaxy{uj)/uj\uj^o is strongly enhanced due to the 
CVC, which is consistent with the analysis in the previ- 
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FIG. 4: Obtained for LSCO by the CVC-FLEX ap- 

proximation. The Drude weight is slightly reduced by the 
CVC. 

ous section. Later, we will discuss its temperature de- 
pendence in more detail. The overall behavior of <Jxy {y>) 
for YBCO is qualitatively similar to that for LSCO. For 
LSCO at r = 0.02, \v&Uxy{,^) takes the maximum value 
at ujxy ~ 0.01, which is about six times larger than ujxx 
for \vi\a(uj). 

The deviation of axy{uj) from the ED-form gives a 
prominent cj-dependence of the optical Hall coefficient 
Rii{u}) — axy{uj)/a'^{uj), which is shown in figlHl For 
LSCO, the w-dependence of R^{uj) given by the RTA, 
axy{uj) / {cr'^ (uj))'^ , is very weak, and its imaginary part 
is tiny. This fact gives the conclusive evidence that 
both cr^y{uj) and (T°(w) follow the ED-form. On the 
other hand, R-ii{uj) given by the conserving approxima- 
tion shows prominent frequency as well as temperature 
dependences. For LSCO at T = 0.02, Imi?H(w) takes 
the maximum value at wrh ~ 0.01, which is similar to 
ujxy for Imaxyiuj) and is six times larger than uJxx for 
Imcr(a;). The relation uJxx S> ujxy,uj-[{}i obtained in the 
present study, which is consistent with experimental ob- 
servation [ia,ll3j cannot be reproduced by the RTA: It 
can be explained only when the back-flow is taken into 
account. Qualitatively similar results are obtained for 
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FIG. 5: Obtained a{u) for YBCO by the CVC-FLEX ap- 
proximation. 



YBCO, although its anomahes are more moderate. The 
observed i?H(^) in YBa2Cu307 at 95K in ref. looks 
similar to the present result for LSCO at T = 0.04 in 
figlHl The model parameters for YBCO used here may 
not really appropriate for a quantitative study. 

In order to elucidate the reason why cra;j,(a;) deviates 
from the ED-form due to the CVC, we analyze Oxyi^^^ 
in the low frequency limit. In the RTA where CVC is 
absent, relations a^y cx and a^y{u})/iijj oc z~'^^^'^ 
are expected. As shown in fig|51 (a), following relations 
are held by the RTA in the present numerical study: 



const., 



ImCT° (wo)/(cr°)^ cx z-^(xT- 



(73) 
(74) 



where luq is a small constant (ujq ^ lO^**). These temper- 
ature dependences is drastically changed due to the CVC, 
as discussed in the previous section. Actually, when the 
CVC's are fully taken into account, we obtain 



/ 2 rp — 0.9 

axy/cr cx a cx J , 
l'm(Txy{uJo)/cr^ cx z^^bocT^ 



-1.7 



(75) 
(76) 
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FIG. 6; Obtained a^y(Lo) for LSCO by the CVC-FLEX ap- 
proximation. The hight of the Drude weight is remarkably 
magnified by the CVC. At T = 0.02, Keaxyi^^) is negative 
for uj > 0.3 due to the CVC, which is consistent with experi- 
ments. 



eq. (|63|l . Thus, both of which are enhanced by the 
CVC's as the temperature decreases. We find that 
lmaxy{i^o) / (J^ cx (axyjcP'Y is approximately realized in 
the present CVC-FLEX approximation. 

In the previous section, relations a cx and 6 = ^™ 
(to = 2 ~ 4) are derived from the analysis of the CVC. 
By eqs.lO-Gnil, relations « T°-^ and to « 3 are 
concluded. To confirm these results more completely, we 
perform another plot shown in fig. |^ (b). 



Oxylaly cx ^2 cx T 
Imcra;j,(a;o)/Imo-°y(LJo) cx ^" cx T 



-0.85 



(77) 
(78) 



where coefficients a and h had been introduced in 



As a result, the relation m « 3 is also derived from 
eas. (|77|) and H78I) . Note that the exponents of T in 
eqs.((7Sl) and ((TTj) are shghtly different because ((t'')^/(t°j^ 
shows a subtle temperature dependence. 

According to eq. lim^^^o ^H('^)/jw = 2z~^{b — 

a)(27kc)~^- As we discussed above, b—a is positive and is 
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FIG. 7: Obtained a^y{Lj) for YBCO by the CVC-FLEX ap- 
proximation. 



enhanced as T decreases. As a result, in nearly AF Fermi 
liquid, icj-linear term of i?H(w) is strongly enhanced by 
the eye, which is consistent with experiments 0, 0, 

mill Hi. 



C. /-sum rule 

The /-sum rule for a^^{uj) gives a rigorous rela- 
tion between the conductivity and the electron density 
liHilli^. It is violated in the RTA because the conser- 
vation laws are not satisfied. On the other hand, /-sum 
rule is automatically satisfied in the conservation approx- 
imation, if all the CVC's given by the Ward identity are 
taken. Thus, /-sum rule is a useful check for the reliabil- 
ity of the numerical study. 

The /-sum rules for a{Lo) and axy{iij) in an anisotropic 
system are given by 



(iwReCT(aj) — ire 



k 



(9fc? 



■7^k, 



(79) 



FIG. 8: Obtained Rii{uj) = (j^y{uj) / {o^y{uj)f by the CVC- 
FLEX approximation. Im_RH(<^) takes large values because 
the ED-form for axy{ijj) is violated due to the CVC. We note 
that the magnitude of in this figure is too large since 

the renormalization to the QP velocity by dSi^/dkx has been 
dropped. The correct value is -Rh""^^ ~ 1 according to ref.^^l- 



/•OC 

/ dujKeaxy{oj) = 0. 
Jo 



(80) 



Equation (|79|l . which can be derived directly from the 
Kubo formula represents the contribution by the 
diamagnetic current. Equation H80() is easily recognized 
form the fact that axy{(^) ^ \^\^'^ a-s \uj\ — > oo, and it is 
analytic in the upper-half plane of the complex w-space. 
In the case of t' = t" = 0, the right hand side of eq.O 
is equal to — (7re^)(ek), where (ek) = X^k ^k^T-k gives the 
kinetic energy. 

The numerical check for the /-sum rule is shown in 
Fig. ^1 Sum{Recr(a;)} and 7r(eJ^) represent the left- and 
right-hand-side of eq.ljTSJ, respectively. Sum{Recr(w)} 
is obtained by performing the numerical w-integration 
form to 100. We see that the /-sum rule ((7^ holds 
well, within the relative error ~ 2%. This results as- 
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FIG. 10: /-sum rules both for a{u!) and axy{uj) are well 
satisfied if the CVC is correctly taken into account. This fact 
assures the reliability of the present numerical study. 



The realization of /-sum rules confirmed in the present 
numerical study is better than expected, despite that all 
the AL-type vertex corrections are dropped. This result 
strongly suggests that the AL terms are insignificant for 
the quantitative study of cr(w) and axy{uj), as they are 
for cr(0) and a^yiO) Ej. 



D. Inverse Hall Angle 



FIG. 9: (a) a^y/a'^y oc r"° *^ at = and Ima^y/lma^y oc 
r~^'^ for small ljq- This result suggests that they are pro- 
portional to 5^ and being the AF correlation length) 



for T = 0.03 - 

Imaxy/u^y oc T' 



' 0.1. (b) o^yla"- oc T" 
'^■^ for small loq. 



at a; = and 



sure the high reliability of the present numerical study 
when Lo is not so large. In general, the Pade approxi- 
mation for larger lo is less reliable because the distance 
form the imaginary axis is large. On the other hand, 
Sum{Recr«TA(^-)| ^-^^^^^^ the RTA (without any CVC) 

is smaller than the correct value, whose relative error is 
more than 12%: This discrepancy is due to the violation 
of the conservation laws in the RTA. 



X 



X 



We also plot / duS^&axy(yi)l / duj\^idijj^y{uj)\ in 

__ Jo "'0 

fig 1 101 where we put X — 100. It should vanish iden- 
tically when X = CO according to the /-sum rule (I80f) in 
the conserving approximation. It becomes less than 0.02 
as shown in fig^l which also suggests the high reliability 
of the present numerical study. This result means that 
the unessential poles of Uxy(y>^ in the upper-half-plane of 
the complex w-plane, which arises from = G^G^ m 
the presence of interaction, are correctly cancelled by the 
vertex corrections. 



IR optical Hall angle Qyl{^) — cr^^y (w)/o-(tj) 
{uj < lOOOcm^^ = 1440K) has been intensively mea- 
sured by Drew et al 0, 0| . They concluded that (I) 
l-m6}i{uj)/L0 is almost independent of w and T, and (II) 
Re0H(w) is also independent of w, while its T-dependence 
is large. In contrast to (I), |Imcr(a;)/a;| monotonously de- 
creases as u! increases, as shown in fig. |31(a). As a result, 
the Hall angle in HTSC follows a simple Drude expression 
IR range (w < 1000cm~i): 



nt 



(81) 



1.5 - 



where and fl^ are w-independent, In contrast, 7ed('-^) 
deduced from the optical conductivity is approximately 
7ed(w) oc maxjtj, ttT} 23]: It is proportional to lo and is 
temperature-independent has a large w-dependence when 
uj > ttT, as recognized in fig. |31 f^Jj is a constant inde- 
pendent of LO and T, whereas it increases as the doping 
decreases. This unexpected behavior of the Hall angle 
puts very severe constraints on theories of HTSC. 

From now on, we show that such anomalous behav- 
iors of 0h(<^) in HTSC's are well understood in terms 
of the Fermi liquid with strong AF fluctuations. The 
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FIG. 11: Obtained tj-dependence of the inverse Hall angle for 
several temperatures. We see that ReS^^ is approximately lj- 
independent for u; < 0.2 while it's temperature dependence is 
large. This result, obtained by taking the CVC into account, 
is consistent with experiments. 



FIG. 12: Obtained temperature dependence of fl^/'y^ = 
9^^{0) and = —ImO'^^ (lj) / lo for several a;'s. Q,^ is almost 
independent of ui and T. 7h is insensitive to uj, whereas it 
strongly depends on T. These behaviors of fl^ and 7^ are 
the most prominent experimental results. 



frequency dependence of back-flow is crucial to repro- 
duce the correct results. Here, we mainly show nu- 
merical results only for LSCO, although similar results 
are also obtained for YBCO. In the present numerical 
study, we derive 0h(w)'s from the analytic continuations 
of Kxy{iuJi)/{K^xii^i) - Kx3:{0)) with uji > 0, which is a 
analytic function on the upper- half complex cj-plane |43j| . 
The value of 0h(w) obtained by this procedure is more 
accurate than dividing axyioj) by (j{lo) after the analytic 
continuations of Kxy{iuji) and K^xii^i) individually. 

Here we discuss the inverse Hall angle by the CVC- 
FLEX approximation for LSCO in more detail. Fig- 
ures ^2 shows the w-dependence of 6^^{uj). We can 
see that Ke9^^{u>) given by the CVC-FLEX approxi- 
mation is almost w-independent for lo < 0.2, which is 
the main experimental finding as explained above. On 
the other hand, Ke{6'^{uj)}^^ by RTA shows sizeable 
w-dependences, which is proportional to 7av(w) given 
in eq.lO. We stress that the effect of the CVC on 



Re^jj (w) is prominent till frequencies much larger than 
7(0). 

Figure ^2 shows the imaginary part of the inverse Hall 
ang le. Im9^\uj) given by the CVC-FLEX approxima- 
tion shows an almost complete w-linear behavior, and 
its gradient stays unchanged against the temperature 
(T = 0.02 ~ 0.1). On the other hand, Im{e^{uj)}~^ 
by the RTA shows a sub-linear behavior with respect to 
u) as shown in fig. 1111 which is inconsistent with experi- 
ments. Its temperature dependence is also inconsistent, 
which will be shown in fig. 1121 Such excess u- and T- 
dependences of {^h(^)}^^ by RTA come from z~-^{uj), 
which will be discussed below. 

Here we further analyze the temperature dependence 
of the inverse Hall angle to make comparison with exper- 
iments. Figure El shows obtained Re6'jj^(ti;) = i}^/2j'^ 
and (— Im0y^(aj)/a;)~"'^ = il^. Apparently, both quanti- 
ties are almost w-independent for lu < 0.2, which confirms 
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the experimental simple Drude expression for the Hall an- 
gle. Here, the value of d in Re^jj^ for w < 0.2 is « 1.6 
in the present study for LSCO (n = 0.9), whereas c? w 1 
by the RTA. Experimentally, d ~ 1 for a; ^ lOOOcm^^ 
in the optimally doped YBa2Cu306+a; {x = 0.93), and 
d ~ 2 in a slightly under-doped compound {x — 0.65). 
Moreover, the value of Re6'j^^(ci;) by CVC-FLEX approx- 
imation is much smaller than that by the RTA, which is 
consistent with experimental observations [TEl [l7| . We 
note that the value of d in the DC-inverse^Hall angle is 
app roximately 2 in under-doped BSCCO 47] and YBCO 
|4q|. It slightly decreases with doping, and d « 1.75 at 
optimum doped systems. 

Figure m also shows that by the CVC-FLEX ap- 
proximation is almost to- and temperature- independent, 
which is consistent with experiments. Contradictory to 
experiments, however, it monotonously increases within 
the RTA. We also stress that the experimental doping de- 
pendence of iljj, which increases as the doping decreases, 
is reproduced well in the present study. According to 
eq. (EZIl, lim^^^o f^n = za"^ {2b - a)'^ . In the RTA where 
a = & = 1, we obtain {£7^}* cx z, which is a increase 
function with T. The inferred T-dependence of H.^ by 
RTA is recognized by the numerical study in figs. Illl and 
[T^ If the CVC's are taken into account, on the other 
hand, the temperature dependence of fl^ will be small 
because fijj ~ za^b~^ is almost constant according to 
eqs.ldl) and In fact, by the CVC-FLEX ap- 

proximation is insensitive to T and uj as shown in fig. 1121 
Experimental observations in HTSC's support the results 
by the CVC-FLEX approximation satisfactorily. 

Here we discuss experimental behavior of the Hall an- 
gle in HTSC's in more detail. In the IR (w = 900 - 
llOOcm"^) measurement the simple Drude form in 
ea. (|81|l is satisfied very well. It is also well recognized 
in YBCO, however, the extrapolation of Im6'jj^ to cj = 
gives a positive intercept, which is recognized as a conse- 
quence of the chain contributions tooxx in YBCO. Cor- 
responding to this fact, reference jig reports that the 
far-IR (w = 20 - 250cm-i) Hall angle in YBaaCuaOr 
deviates from ea. (|FT|l . One possible origin of this de- 
viation other than the chain contribution would be the 
emergence of the pseudo-gap. In fact, DC transport coef- 
ficients show various anomalous behaviors in the pseudo- 
gap region. They are well reproduced theoretically in 
terms of the AF-I-SC fluctuation theory if one take the 
CVC into account ;14 j . It is an important future problem 
to extend the scope of the present study to the pseudo- 
gap region, using the AF-I-SC fluctuation theory. 

In summary, experimentally observed simple Drude 
form of 6'jj^(ti>) in eq.JHIJ is satisfactorily well repro- 
duced by the CVC-FLEX approximation, using the pa- 
rameters for LSCO. Both 7^ and fJJj are constant for 
u; < 0.3, while the former is strongly temperature depen- 
dent. Similar results are obtained even if one use pa- 
rameters for YBCO, as shown in fig. ^| We note that 
= 1.0 in the upper panel of Figll2l corresponds to 
5000 [Tesla/radian], and OJj = 0.2 in the lower panel 
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FIG. 13: Inverse Hall angle for YBCO obtained by the CVC- 
FLEX approximation. 

corresponds to 0.2 [1/cm Tesla], approximately. These 
values seem to be well consistent with experiments [T^ . 



E. Hall Angle 

We also discuss 0h(w) given by the CVC-FLEX ap- 
proximation, and make comparison with experiments. 
Figured shows V^cQy^W) for LSCO. 6*° (w) by RTA is al- 
most temperature independent for w > 0.15. In contrast, 
6'h(w) by the CVC-FLEX approximation is T-dependent 
till much larger uj due to the w-dependence of the CVC. 
Especially, it increases with T for lo > 0.15. which is 
consistent with experimental observation [isl Il7l | . This 
change of dReOn/dT for larger a; is a natural conse- 
quence of the Lorentzian form of Re0H('j-'), where 7h is 
cj-independent and is an increase function of tempera- 
ture. In contrast, Recr(a-') deviates from the Lorentzian 
due to the a;-dependence of ^■ET>{^)^ which is consistent 
with experiments jlSj, iTQ] . 

The origin of the Lorentzian form of YyeO'ii{Ljj) is as- 
cribed to the almost perfect cancellation of w-dependence 
of 7ed(w) and that of the CVC: Rc6'H(ti^) for w » 7 
will be enhances by the former effect because 7ed(w) is 
a increase function of w, whereas it will be suppressed 
by the the latter, because the back-flow will be less im- 
portant for larger w. It is a nontrivial future problem 
why these two effects cancel out almost completely, which 
results in the observed Lorentzian form of Re0H(w) for 
Lo < lOOOcm-1. 

FigurelT^lshows the temperature dependences of 0h(^^) 
for several w's. The obtained results for oj — 0.2 or 0.4 
look similar to the experimental observations for YBCO 
with X — 0.93 (optimum) or x = 0.65 (slightly under- 
doped) for UJ ~ lOOOcm"^ [l^. On the other hand, the 
result for lo — 0.04 resembles the observation for heav- 
ily under-doped non-superconducting sample (x = 0.4). 
We guess from this fact that the electronic states in heav- 
ily under-doped systems are qualitatively reproduced, al- 
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FIG. 14: Obtained Li;-dependence of the Hall angle for several 
temperatures. 
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FIG. 15: Obtained temperature dependence of the inverse 
Hall angle for several cj's. Characteristic behaviors of the 
experimental 0h are well reproduced. 



though the experimental value of 7^ is much larger. 



F. Predictions for Electron-Doped Systems 

DC transport phenomena under magnetic field in 
electron-doped systems (e.g., NCCO) also shows striking 
NFL behaviors which originate from the CVC |13, ilJ] ■ 
Surprisingly, both _Rh and S* in a under-doped NCCO is 
negative, and its absolute value increase as T decreases. 
Their behavior looks approximately symmetrical to those 
in hole-doped systems. Contrary to these experimen- 
tal facts, the RTA predicts the positive Hall coefficient 
because it has a hole-like FS whose shape is similar to 
YBCO. This discrepancy is naturally solved if one take 
the CVC into account, since dO'^ /dk\\ becomes positive 
around the cold spot of NCCO whose location is different 
from that of YBCO; see figll](a) 1^. 

Quite recently, optical Hall conductivity in electron- 
doped systems has been observed by Zimmers et al 
|2l| . Here, we analyze the w-dependences of a^^{uj) in 
electron-doped systems based on the conserving approx- 
imation. Figure [TCI shows cr(Lj), Uxyioj) and i?H(t^) ob- 
tained by the CVC-FLEX approximation. Both (Jxy{uj) 
and Ryl{uj) for NCCO are similar to those for LSCO 
given in Figs. and |H1 except their signs. We stress 
that Imi?H(w) is as large as Rei?H(i^) for finite w, which 
means that the simple ED-form of axy{ijj) is violated. We 
predict that the signs of Rei?H('^) and Recrxi/(w) change 
from negative to positive with lo. Thus, the CVC in 
NCCO plays important roles. In future, measurements 
of axyioj) in NCCO are highly anticipated. 

We found that an accurate numerical calculation (Fade 
approximation) for NCCO is much difficult than that for 
LSCO and YBCO. By this reason, we could not obtain 
reliable results for 0.04 >T> 0.08. It is a future impor- 
tant problem to improve the stability of the Fade approx- 
imation in case of NCCO. 



V. SUMMARY AND FUTURE PROBLEMS 

In the present work, we have calculated the optical 
conductivities a-{i^) and axy{i^) for HTSC's by the CVC- 
FLEX approximation. Experimentally observed anoma- 
lous behaviors for (t(w), axy{oj), and i?H(w) are 
well reproduced for enough wide range of frequencies 
and temperatures, without assuming any fitting parame- 
ters "22]. Especially, (I) a{u;) given by the CVC-FLEX 
approximation follows the ED-form shown in eq. H7U|) 
with the relaxation time in eq. H72|l . whereas a^yiuj) 
strongly deviates from the ED-form, eq. (|71|) . because 
the w-dependence is much exaggerated due to the CVC in 
nearly AF Fermi liquids. By this reason, (II) ImR}i{Lu) ~ 
Rei?H(w) is realized even when for w ^ 7, as shown in 
Fig. |S1 Moreover, (HI) 0h(w) follows a simple Drude 
form given in ea. (|81|) for uj < 0.2, as shown in Figs. ^2 
[T^ and ^1 They are consistent with the characteris- 
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FIG. 16: Obtained a{uj), a^y{uj) and Rh{uj) for NCCO by 
the CVC-FLEX approximation. 



tic experimental results for HTSC's reported by Drew et 
al. 15, 17, 19J. These anomalous AC transport phe- 
nomena cannot be reproduced by previous theoretical 
works based on the RTA, even if one assume extremely 
anisotropic Tk. Instead, they are naturally explained 
by taking the CVC into account in accordance with the 
Ward identity. 

In the present study, we have pointed out the impor- 
tant role of the back-flow in the optical conductivities for 
the first time. The enhancement of ImtT2;y(a;)/iix'|i^=o due 



to the CVC is not same as the enhancement of axy{0); 
The former is more prominent than the latter as ex- 
plained in eas. (|75|l - l|78|l . This fact leads to the break- 
down of the extended Drudc-form at very low frequencies. 
The back-flow decreases monotonically with uj as one ap- 
proaches the collisionless region (w >, 7). This fact gives 
an approximate Drude-form of the Hall angle in ea. H81|) 
for oj < 0.2, nonetheless of the fact that cr(w) deviates 
from a simple Drude-form (instead it follows a ED-form) 
due to the cj-dependence of ^{oj) . Note that the back-flow 
in the collisionless region is given by the real part of 
PgL I23, |3^ • This is an important future problem for us 
to find a simple physical explanation for this numerical 
result. 

We stress that both AC and DC anomalous transport 
phenomena in HTSC's are explained in a unified way 
based on the Fermi liquid theory, if one take the CVC to 
satisfy the conservation laws. As for the DC Hall coef- 
ficient, one frequently attribute the enhancement of Rn 
to the small area of the cold spot (Fermi arc) observed 
by ARPES in under-doped compounds. However, this 
idea contradicts the fact that the i?H decreases in the 
pseudo-gap region while the Fermi arc shrinks further. 
In the same way, anomalous behaviors of Ap/p, S and 

in the pseudo-gap region cannot be understood within 
the scheme of the RTA. Such contradictions are satisfac- 
torily solved by the CVC-FLEX approximation, by tak- 
ing the superconducting fiuctuations induced by the AF 
fiuctuations Ti"] . We stress that the natural extension of 
this DC transport theory to AC transport phenomena, 
with taking the same diagrams for the CVC, succeeds 
in explaining the optical Hall effect observed in HTSC's. 
This fact means that the qualitative dynamical electronic 
properties of HTSC's, from the over-doped to the slightly 
under-doped systems, are well understood in terms of the 
Fermi liquid theory with strong AF fluctuations. 

There remain many important issues for the Future 
study. For example, one can study various AC-transport 
coefficients other than based on the CVC-FLEX ap- 
proximation, using the similar method developed in the 
present study. Study of the role of the CVC at finite 
frequencies for Ap/p, S and u would be very interesting, 
although experimental observation would be difficult at 
the present stage. In addition, we are planning to study 
the optical conductivities in the pseudo-gap region based 
on the FLEX+T-matrix approximation, which ascribes 
the pseudo-gap phenomena in HTSC's to the strong su- 
perconducting fiuctuations hj. As we mentioned in the 
previous section, reference Tq reports that the far-IR 
(cj = 20 - 250cm~i) HaU angle in YBa2Cu307 devi- 
ates from the Drude-form in ea. H81|) . although it is well 
satisfied for uj ~ lOOOcm"^. We would like to find out 
whether such an anomaly in far-IR Hall angle could be 
understood as the pseudo-gap effect using the AF+SC 
ffuctuation theory. 



20 



Acknowledgments 

The author is grateful to H.D. Drew and A. Zmimers 
for fruitful discussions. 



APPENDIX A: PHYSICAL MEANING OF THE 
BACK-FLOW 

Throughout the present work, we have stressed the im- 
portance of the back-flow for both DC and AC transport 
phenomena. Here, we would like to depict the physi- 
cal aspect of the back-flow in nearly AF Fermi liquid 
based on the phenomenological Landau-Fermi liquid the- 
ory According to Landau, the energy of the quasi- 
particle eko- is expressed as 

Ekcr = CkCT + ^ /kcr,k'cr"5?lkV + 0(((5n)^), (Al) 



where eto- = k^/2m*, nta- is the distribution function 
of quasiparticles, /ko-.k'o-' is the Landau function, and 
Sn^f^a = "-ko- — ^ka^- Equation (|A1I) means that the en- 
ergy of quasiparticles are changed when the quasiparticle 
excitation exists. By this reason, once we add a quasipar- 
ticle at k outside the FS, the Fermi sphere is deformed 
to minimize the total energy unless /k(T,k'cr' =0. As a 
result, the Fermi sphere has a finite momentum, which 
is the physical meaning of the back- flow. Thus, the ex- 
istence of the back-flow is assured by the most essential 
relation in the Fermi liquid, ea. l|Aip . Apparently, the 
back-flow would be indispensable in strongly correlated 
Fermi liquids, like in HTSC. 

The importance of the back-flow has been understood 
very well in a spherical system, where /k<j,k'<T' can be ex- 
panded by Legendre polynomials, P;(k- k'). Landau first 
studied the back- flow in the coUisionless region w ^ 7, 
where the lifetime of an quasiparticles is longer than the 
period of the outer field. Based on the Kubo formula, 
Yamada and Yosida analyzed the opposite region cj ^ 7 
in order to study the role of the back-flow on the DC con- 
ductivity. They rigorously proved that the conductivity 
diverges even at finite temperatures if no Umklapp scat- 
tering process exists. In contrast, the RTA always gives 
finite conductivity cr oc 7^^ at T 7^ even in the absence 
of the Umklapp process, reflecting the the violation of 
the momentum conservation laws. 

In contrast, importance of the back-flow in anisotropic 
systems with strong correlations has not been recognized 
until recently. As explained in §111, we found that total 
current Jk becomes quite different from the quasiparti- 
cle velocity Vk due to the CVC when the AF fluctuations 
with q « Q are strong J^] . This is the origin of various 
anomalous transport phenomena in HTSC's. This unex- 
pected behavior of Jk comes from the fact that (k—k') 
in the Bethe-Salpeter equation (|50|l . which corresponds 
to the Landau function for w <C 7, takes large values 
only for k — k' w Q. In this case, according to ea. ljAljl . a 



quasiparticle added at k' strongly modifies ek only when 
k « k' -I- Q, which makes the induced current (back-flow) 
proportional to Vk. The induced current is not parallel to 
the source velocity Vk' , in contrast to the case of spher- 
ical systems. The schematic behavior of Jk in HTSC's 
is shown in flg^ (b) . Jk at the hot spot takes enhanced 
values because ak <, 1 in ea. H52|) . which is interpreted as 
the "resonance" between Vk and Vk'. 

In the present paper, we studied the optical conduc- 
tivity and Hall conductivity by taking the w-dependence 
of the CVC into account appropriately, which has not 
been performed in previous studies. We flnd that axy{ijj) 
shows a striking w-dependence when the AF fluctua- 
tions are strong, which cannot be expressed by an ED- 
form. Such a non-Fermi liquid-like behavior comes from 
the prominent w-dependence of the CVC, which was 
detected in the present study for the first time. As 
shown in §111, the total current at finite oj is given by 
Jk(w) — 3k + *^Jfe^'' + 0{uP'), where J^"'^'' is real and its 
k-dependence is much larger than the first term. This 
strong w-dependence of the total current gives rich vari- 
ety of spectrum in optical conductivities. 



APPENDIX B: COMMENTS ON PREVIOUS 
THEORETICAL STUDIES 

Anomalous DC transport phenomena in HTSC's, as 
represented by the enhancement of the Hall coefficient, 
have been frequently ascribed to the reduction of the ef- 
fective carrier number within the RTA. For example, Ref. 
[35^1 proposed the highly anisotropic Tk model based on 
a spin fiuctuation theory; Th oc for hot electrons 

whose density is Tih, Tc oc T^^ for cold electrons whose 
density is nc{— n — n^). They assume that ric <C rih and 
Tc/Th ~ 100 at lower temperatures. Their model can- 
not give a comprehensive explanation for anomalous DC 
transport phenomena in HTSC's, while the CVC-FLEX 
approximation can give it. 

Here, we examine the optical conductivities within the 
RTA based on a simplified anisotropic 7k model as fol- 
lows: 



a{uj) cx 



27c ~ ioj 27h — ioj 
ric , nil 



(27c - iuj)^ (27h - iiuy 



(Bl) 
(B2) 



where rh,c = l/27h,c- The Hall coefficient is highly en- 
hanced in proportion to l/enc(oc T~^) when 7ic/7c ^ 
"h/7h- 

In the case of ric <SC rih and 7c <C 7h, frequencies lJxx, 
LOxy and o^RH which give the maximum Imcr(Li;), Im(Tj;j,(a;) 
and Imi?H(w) respectively, are given by 



^xy 
Wrh 



27c, 

(2/x/3)7c = 1.167c 
2(7cnh + 7h"'c) 



nil 



(B3) 
(B4) 

(B5) 



21 



where cjrh is approximately given by 

LJRH ~ 1.16nc7h/nh > 7c for nc/7c > rih/7h, 
~ 2.3l7c for nc/7c = nh/7h, 

~ 1.167c for nc/7c < nii/7h. 

Because nc/7c(oc T"'^) will be larger than nh/7h(oc T~^) 
at lower temperatures, the relation o^rh ^ i-^Ka; is ex- 
pected in this model. This result is inconsistent with 
experimental fact wrh ^ i-^xx, as mentioned in §1. 

In a similar way, d-density wave (DDW) model j50l | 
have been proposed to explain the enhancement of of the 
Hall coefficient: Rn increases below the d-density wave 
transition temperature, inversely proportional to the area 
of the "Fermi arc" . However, it is hopeless to reproduce 
the characteristic experimental behavior of Rb_{uj) in this 
model. 

We also comment on the 2D Luttinger liquid model 
with two kinds of relaxation times (rtr cx: T~^, th oc T~^) 



proposed by Anderson [sj]. In this model, DC con- 
ductivities are given by cr oc and axy oc (rtrTn)"^, 
respectively. Their natural extensions to the optical 
conductivities are given as 17(0;) oc (r^~^ — iuj)~^ and 
axy{uj) oc (rtr/TH)cr^(w) j^. This result directly means 
that Imi?H('-^) = 0, which apparently contradicts ex- 
periments. In addition, the Hall angle in this model is 
6'h('^) oc (Ttr/Tn) • (T(.~^ — iu})~^] the temperature depen- 
dences of coefficients are different from the experimental 
ones given in ea. (|81|l . Note that another functional form 
of (7^y{uj) which predict finite Imi?H('^) was proposed in 
ref. (I9j, although its theoretical verification is uncertain. 

In summary, a comprehensive understanding for the 
optical Hall coefficient in HTSC's cannot be obtained by 
previous theoretical works based on the RTA, or by the 
2D Luttinger liquid theory. The transport theory based 
on the Fermi liquid theory presented in the present pa- 
per, where the CVC is correctly taken into account, can 
explain various experimental anomalies at the same time. 
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